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Abstract. Magnetoplasmadynamic (MPD) thrusters are capable of accelerating quasi-neutral plasmas to high exhaust 
velocities using Megawatts (MW) of electric power. These characteristics make such devices worthy of consideration 
for demanding, far-term missions such as the human exploration of Mars or beyond. Assessment of MPD thrusters at 
the system and mission level is often difficult due to their status as ongoing experimental research topics rather than 
developed thrusters. However, in order to assess MPD thrusters’ utility in later missions, some adequate 
characterization of performance, or more exactly, projected performance, and system level definition are required for 
use in analyses. The most recent physical models of self-field MPD thrusters have been examined, assessed, and 
reconfigured for use by systems and mission analysts. The physical models allow for rational projections of thruster 
performance based on physical parameters that can be measured in the laboratory. The models and their implications 
for the design of future MPD thrusters are presented. 


INTRODUCTION 

Magnetoplasmadynamic (MPD) thrusters have long held the promise of high exhaust velocity (u e ) at MW power 
levels (Sovey, 1991). The combination of high u e and high power in a compact device is especially beneficial for 
demanding missions such as the human exploration of other planets, which will require lightweight, high power 
density propulsion to be feasible. At present, MPD thruster research is still at a fundamental rather than 
developmental level. Research goals include increasing thruster efficiency and operating lifetime. It may be safely 
said that current experimental MPD thrusters do not represent the refined and optimized thrusters that might be used 
in missions planned a generation from now, so curve fits to present data will not accurately characterize the MPD 
thruster of the future. In order to accurately assess their promise, a method of characterizing future thruster designs 
based on current knowledge of their physical behavior is needed. These characterizations are needed by system and 
mission analysts to assess MPD thrusters’ benefits to later missions. 

A secondary consideration for practical assessment of MPD capabilities is the translation of experimental data into a 
form suitable for system integration and parametric studies. In the laboratory, an experiment is generally operated 
with the thruster current and mass flow rate as the control variables, and thruster voltage and thrust as the 
measurements. From these variables, u e , power (P e ), and thrust efficiency (r|) can be calculated. In contrast, a 
mission analyst uses u e , P e , and r| to determine the trajectory and performance of a low-thrust spacecraft. A system 
analyst needs to derive the current, voltage, and mass of the MPD thruster from the mission u c and P e in order to 
then design the power system and spacecraft. Thus, MPD thruster models developed to explain thruster behavior in 
the laboratory must be converted to a variable space of use to system designers. 

MPD experiments and theories describe MPD thrusters using variables that are relevant to experiment such as mass 
flow rate ( m ) and current (J): 



( 1 ) 


V = V(m,J ) ( measured ) 

T = T(m,J ) ( measured ) 
u e = u e (m,T) (calculated) 
r\ = i](m,J,V) ( calculated ) 

with the thruster geometry and plasma physics included in the relations for V and T. The system model reduces the 
number of free variables to only 2, from which all other parameters are derived: 

V=V(u e ,P e ) 

J=J{u e ,P e ) 

m= m(u e ,P e ) ( 2 ) 

?7= r] (u e ,P e ) 

T = mu e 

It is important to note that in order to derive the system equations from the experimental ones, the models must be 
mathematically structured to allow analytical conversion from one “frame” to another. This is an artificial 
requirement imposed by the need for simple models to be incorporated into large, complex system or trajectory 
codes. 

The analysis described herein consists of several parts. First, existing self-field MPD thruster models were derived 
using recent models for thrust production in conjunction with a voltage model that includes various known energy 
channels available in an MPD thruster discharge. Second, the models considered were converted using 
Mathematica from the experimental variable frame to the systematic frame. Any “free” parameters in the models 
were identified. Third, the models were compared to experimental data for two hydrogen MPD thrusters. From this 
data, free parameters were determined and used to assess the model over a range of operating conditions. 


MPD Thruster Description 

The Magnetoplasmadynamic thruster, or MPDT, is a coaxial Lorentz force accelerator. A schematic of self-field 
MPD operation is shown in figure 1 . The thruster consists of an outer annular anode, a central cylindrical cathode, 
and an insulated backplate. A discharge is created between the electrodes, with current flowing from anode to 
cathode. The self-induced magnetic field from this current is primarily azimuthal, resulting in a j r X B„ Lorentz 
force acting on the quasi-neutral plasma in the axial direction. A thruster using only the self-generated field is 
referred to as a “self-field” MPDT. An additional, externally generated axial magnetic field may also be imposed on 
the thruster in an effort to increase performance; this is referred to as an “applied-field” MPDT. The anode is 
typically constructed of a high temperature refractory metal such as tungsten or molybdenum. The cathode is 
tungsten impregnated with thorium or barium oxides to reduce the work function for electron emission. The 
insulating backplate is boron nitride. Propellant is injected at the rear of the thruster, or possibly along the sides, and 
ionized by the current flowing through the plasma in the thruster. Because the self-field thruster performance 
depends on the current used in the thruster, high (kiloamperes) currents are required for high specific impulse. 
Because the propellant must be ionized in order for the thruster to operate, propellants with low first ionization 
energy and high second ionization energy are predicted to give the highest efficiency. Propellants of choice for 
MPD thrusters are hydrogen and lithium, because of their low atomic mass and relatively low ionization energies 
(Ageyev, 1993 and Lapointe, 2002). Lithium’s very low first ionization energy makes it particularly advantageous 
at lower powers (0.5 - 5 MW) and moderate exhaust velocities (20-70 km/s). Hydrogen is more amenable to 
multimegawatt and high specific impulse operations, is easier to supply to the thruster as a gas, and is less likely to 
condense on spacecraft surfaces. 
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FIGURE 1. 


Diagram of MPD Thruster Processes. The radial current and azimuthal magnetic field produce thrust. 


SELF-FIELD MPD THRUSTER MODEL 


Although conceptually and often geometrically simple, MPD thrusters have been found to be rich in the complexity 
of the plasma flows they generate. There are several MPDT thrust models in the literature; of these, the most recent 
self-field models were examined for this system model: 1.) the modified Macker model, (Jahn, 1968, Choueiri, 
1987) 2.) the Tikhonov model (Tikhonov, 1993), and 3.) an empirical model by Choueiri (Choueiri, 1998b). Each 
of these models have been used to model hydrogen MPD thruster data. For the model comparisons, each thrust 
model has been coupled to a single common voltage model. The relative strengths and weakness of each model will 
be evaluated, and improvements will be suggested. 


The voltage model used with all three thrust models considers various energy partitions possible in the plasma.: 
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The first term represents the directed kinetic energy voltage due to the acceleration of the plasma; a divergence term 
(5) was added to represent the thrust losses due to the divergence of the plasma flow, which reduces effective thrust. 
The power required to ionize the plasma is modeled as an ionization energy per kilogram (e;) for the propellant. 
Depending upon the MPD thrust model assumed, the ionization fraction may vary with the exhaust velocity or 
current. Specifically, some models have the ionization fraction linearly proportional to the current at low currents; 
this variation is signified by the gj term. A heuristic ohmic heating term (f J/a) is included, with the geometric factor 
“f ’ to account for varying thruster geometries. In later comparisons, a o value of one-tenth Spitzer conductivity was 
used to account for anomalous resistivity observed in MPD thrusters (Choueiri, 1993). The geometric term, f, was 
approximated as ln(r a /r c )/(2:tL) to give some variation in ohmic losses depending on thruster dimensions. Finally, 
electrode and sheath potentials representing power deposited to the cathode and anode are represented simplistically 
by a constant “fall” voltage, Vf. This term is the only “free” parameter in the models (apart from ionization loss, to 
be discussed later) consists of the sum of the plasma sheath voltages at both anode and cathode, and was observed to 
be thruster dependant. The term ranges from 50 - 150 V and is empirically determined to best fit the data. 
However, experimental tests have demonstrated (Myers, 1992) that there is a physical basis for this fall voltage 
term, and that each thruster has its own characteristic fall voltage. 



Modified Macker Model 


The Macker formula is the most basic description of thrust generation by an MPD thruster. The relation calculates 
the J X B Lorentz force for coaxial cylindrical electrodes, assuming azimuthal uniformity, purely radial current 
except at the cathode tip, and a simplified axial distribution current along the cathode. The model takes the form 


T=bJ 2 


( 4 ) 


The coefficient “b” is related to the thruster electrode geometry: 


b = 


V o 
4jt 



\ 

+ a 
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( 5 ) 


where r a is the thruster anode radius, r c is the central cathode radius, and a is a term accounting for thrust 
contributions from the cathode tip (Jahn, 1968). This model was initially used by researchers to describe MPDT 
thrust over the full range of currents and mass flow rates; however, subsequent thrust measurements showed that the 
thrust departed from this model at low values of current and exhaust velocity (Gilland, 1987). A later theory by 
Choueiri (Choueiri, 1987) divided MPDT thrust behavior into super- and sub-Alfven critical speed operating 
regimes, with different thrust - current relationships in each regime. The Alfven critical speed (Ua) is the velocity 
where the kinetic energy of ions in the propellant is equal to the first ionization energy of the neutral atom: 


U A 



( 6 ) 


According to this model’s interpretation, the propellant is only partially ionized in the sub-Alfven regime. The 
model assumes that there is a linear relation between current and ionization fraction in the plasma up to the current 
and mass flow rate at which the plasma is fully ionized, Jg: 


J fi - 



( 7 ) 


This allowance for different thrust profiles at super- and sub-Alfven critical speed necessitates a continuous, 
piecewise function for thrust: 


T(c) = 


m Ua 

b J 


Jfi 

2 


u c < Ua 
u e > Ua 


( 8 ) 


The modified Macker thrust model developed for this paper takes the traditional form at super- Alfven critical speeds 
but uses the partially ionized model below. Because of the piecewise nature of the thrust model, the voltage model 
was defined for the two different regimes to accommodate the partially ionized plasma at currents below Jg: 
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Tikhonov Model 


The Tikhonov (1993) thrust model is very similar to the modified Macker model discussed above. The primary 
difference between the two models is that the Tikhonov thrust model attempts to operate with only one thrust 
equation over the full range of current/exhaust velocity: 


T = 


I b + 


(J 2 Ini) 2 


( 10 ) 


The parameter y is a geometry and propellant dependent constant. The y term in the thrust model decreases in 
magnitude as the discharge current increases, and so there would be a major divergence between the Tikhonov and 
the modified Macker models only at very small J or equivalently, u e . In fact, the y term begins to have an impact on 
net thrust for H 2 propellant only very far below sub-Alfven critical speeds. 


Empirical Model 

The empirical model takes an approach similar to the Tikhonov model, but relies on phenomenological observations 
to match the experimental data, rather than fundamental theories (Choueiri, 1998a): 


T = 




+ 


In 



( 11 ) 


The dependence on Ua is maintained by the term § = J/Jg. The increase in T/J 2 at low currents is modeled by the 5 4 
term, v is an empirical constant that appears to be dependent on propellant species and mass flow rate. It should be 
noted that Choueiri found the original model did not accurately reproduce experimental results with hydrogen, 
possibly because of the variability in defining the Alfven critical speed and therefore the § parameter for a molecular 
propellant. 


RESULTS 

Since these models take a systems analysis point of view, each model had to meet two criteria. First, that the 
equation set could be analytically transformed into system variables. Second, that the resulting model could match 
experimental efficiencies and voltage vs. current data. Using Mathematica, the three models were converted to the 
systems frame and compared to data from the Princeton Benchmark thruster (Choueiri, 1998b) and the Osaka MYII 
thruster (Tahara,1985), varying the free parameter of Vf to best fit the experimental curves. 

Upon mathematical evaluation, the Tikhonov model was not useful because it failed to produce real (non-imaginary) 
results at sub-Alfven critical speed if the fall voltage was kept constant over the whole range of specific impulses. 
Assuming that the fall voltage parameter was allowed to change, the model was able to produce real data for sub- 
Alfven critical speeds, but was unable to accurately track the Princeton benchmark data that was used as 
comparison. Thus, to use this model effectively for systems analysis, one requirement would be a new voltage 
model that could account for a variable fall voltage. In addition, the overall effect of the y term was to scale the 
thrust with J“ 2 only at low u e (< 10 km/s), an effect that is difficult to observe in any experimental data due to the 
difficulty of operating at very low discharge currents. From the systems analysis perspective, this effect was 
determined to be insignificant; few missions can be envisioned which require operation at the u e ’s at which the 
additional term becomes important. For these reasons, as well as the fact that the Tikhonov expression so closely 
tracked the Macker model at super-Alfven speeds, and because of the extremely long compiling times needed to 
reliably run the model, this model was not selected for further testing. 

The empirical thrust model in conjunction with the assumed form of the voltage equation was found to be non- 
analytic in terms of converting from an experimental to a system level formulation. In an attempt to approximate 
the model to allow for the experiment-to-system conversion, a first order Taylor expansion of the original thrust 
model was tried. This approximation was in good agreement with data presented by Choueiri (Choueiri, 1998b) and 



reproduced the relevant features of the model. The empirical mode proved to be virtually identical to the modified 
Macker model over a wide range of power and u 0 ; primary differences occurred at u e less than 10 krn/s, which is 
generally below the range of interest to mission designers. The model’s inaccuracy in reproducing hydrogen data 
resulted in its rejection as a basis for a system model. 

The Macker model was convertible to a system frame, although the resulting models are piecewise continuous and 
must be defined in two regimes. The representative equations from the modified Macker model for J(u e ), 
V(u c ),T(u e ), m(u e ) and r|(u e ) are shown here for the two operating regimes described previously. 

For u e < Ua, the model equations are given with the annotation “1” to denote sub-Alfven critical speed operation: 
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2 b Ua cm 2 b Pe ei a + uil f Pe + (Vf 2 + b Pe Ua (1 + d)) cr) - X 
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where 


b u e Ua cr((2 b Pe ei o + u e ( 2 f Pe + (Vf 2 + b Pe Ua (1 + <5)) a] - x 
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X = -yjiu Vf -Jo 4 b Pe ei o + u e (4 f Pe + (Vf 2 + 2 b Pe Ua (1 + 6)) 
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Similarly for u e >UA, using the annotation “h”: 


-u Vf a + Ju e -\lo t/ 2 b u e ’ Pe (1 + 5) o + 4 b Pe ei a + u e (4 f Pe + Vf 2 o) 


2 u e f + b u e "(l +8) o + 2b ei o 


Jh( u e ) = 


(17) 




Vh( u e ) = 
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( 20 ) 


where 


ip = 


V^VfV^ -^2 b u e 2 Pe (1 +8) 


+ 4 b Pe ei o + u e (4 f Pe + Vf 2 o) 


( 21 ) 


In each case the fall voltage for each model was adjusted so that mass flow, current, voltage, and efficiency all fit 
best at a reference point; however, once the fall voltage was determined for each model and thruster, it was left 
constant over the full course of operation. P e and u c were input directly from experimental data. One outstanding 
issue in modeling thrust was the value of a to be used in the thrust coefficient definition. Common values range 
from 0.25 to 0.75 (Jahn, 1968). In the course of this study, values of 0.0 to 0.75 were considered; a value of a=0.75 
was found to give the best agreement with data. Table 1 indicates the relevant parameters used in the simulations. 
The two thrusters considered are shown in Figure 2. 


TABLE 1. Relevant Modeling Parameters for MYI1 and Princeton Benchmark (PBT) Thrusters Considered in System 
Model. 


Thruster 

r a /r c 

Prop 

Macker V f 

MY II 

3.8 

h 2 

80 

PBT 

5.4 

h 2 

150 


The modified Macker model matched PBT experimental data well for a wide range of specific impulses. Figure 3 
shows a comparison of the modified Macker thrust model with experimental data from the two thrusters considered 
in this study. The “effective b” on the abscissa is simply thrust divided by the square of the current. In the case of 




FIGURE 2. Schematics of a.) MYII Thruster (Tahara, 1985) and b.) Princeton Benchmark Thruster (Choueiri, 1998b) 

(Drawings are not to scale). 
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FIGURE 3. Comparison of Modified Macker Model with Thruster Data at both Low and High Power Levels. 


the model, this value is by definition a constant above Jf„ and follows a J" 1 behavior below Jfi. The data for the two 
thrusters show similar but not identical increase with low J. As has been observed by several researchers (Choueiri, 
1987; Fischer, 1985), the exact value of the Alfven critical speed is difficult to specify for a molecular propellant, so 
that the point at which an actual thruster transitions from fully to partially ionized may not be clearly marked. The 
thrust behavior near Jf, may therefore not be strictly proportional to J" 1 . However, the model and thruster behavior 
appear to follow each other qualitatively in the partially ionized regime. The difference in thrust versus current 
stems from the failure of the full model (thrust, voltage, current, mass flow rate) to simultaneously match all of the 
experimental parameters. 

A direct comparison of model and data voltage-current characteristics for the PBT is shown in Figure 4, and of 
efficiency-Ue in Figure 5. The change in slope of the curves is the point at which J=J ri . A difficulty in comparing the 
model with PBT hydrogen data is the large voltage fluctuations observed in the experiment. The values used in 
Figures 4 and 5 are the mean values. The model follows the current-voltage behavior of the Princeton thruster 
reasonably well, particularly at high current levels. The model appears to over-predict thruster efficiency by 
approximately 5% (i.e., the model efficiency in percent is 5% greater than the measured values). Figure 5 represents 
the data that would be used in calculating optimal trajectories; for scoping and system assessments, the model 
appears to duplicate MPD behavior reasonably well. 

Calculations for the MYII thruster using hydrogen propellant yielded less accurate results. Thruster and model 



FIGURE 4. Hi PBT Terminal Voltage Characteristics and 
Model Results for Two Mass Flow Rates. 


FIGURE 5. Hi PBT Mission Performance Characteristics 
and Model Results for Two Mass Flow Rates. 
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FIGURE 6. MYII Terminal Characteristics (Tahara 1985) FIGURE 7. MYII Mission Performance Model Shows 
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voltage characteristics are shown in Figure 6. The model diverges from experimental behavior by as much as 25% 
at low currents and powers, but closely agrees near and above Jg, which would be 19 kA for this thruster and mass 
flow rate. Because the MYII has a smaller discharge chamber than the PBT and was run at higher mass flow rates, 
additional ionization, excitation, and recombination processes may be occurring in the thruster, particularly at low 
u e . Specifically, recovery of some of the ionization energy from the plasma flow would reduce the amount of frozen 
flow energy convected with the plasma and result in a lower voltage for a given u e . This effect has also been 
observed in MACH2 calculations simulating the same thruster (Mikellides, 2001). In contrast, the modeled mission 
performance shown in Figure 7 agrees quite well with the data and could also be used in a mission analysis of this 
type of MPD thruster. Flowever, it is clear that the processes occurring below Ua in the MYII thruster are more 
complex than assumed in the current model. 


CONCLUSIONS 

The model described above offers some initial descriptions of MPD thrusters on a system level, and provides some 
insight into some of the physical processes occurring in the MPD plasma. The first goal of this study, to generate a 
MPD thruster performance model suitable for use in system or mission analysis, has been achieved. More important 
than the accuracy of the model in reproducing experimental results is the fact that the model can be modified in 
terms of scientifically relevant parameters, such as fall voltage, ionization energy of the propellant, and thruster 
dimensions. This will allow the user to assess the benefit of fundamental improvements to MPD thruster behavior, 
such as reducing the fall voltage, to project possible performance improvements over a range of powers and exhaust 
velocities. 

In comparing and deriving the parameters needed to match the model to experimental results, some insight has been 
gained into the validity of the assumptions. The assumed form of the voltage model used in the set of equations 
imposes the greatest uncertainty in deriving a unified model of the MPD discharge. Of the assumed constant 
parameters in the model - £;, b, Vf— two arise from the voltage equation. The relative fidelity of the model to PBT 
data compared to that of the MYII data implies a marked difference in the behavior of these two devices over a 
range of power. The simple assumptions of constant ionization energy, constant fall voltage, and constant thrust 
coefficient followed PBT behavior relatively well. Some or all of these assumptions break down when modeling the 
MYII. 

Some additional aspects of MPD performance have not yet been addressed with this model. Not yet assessed in this 
model is the effect of propellant species. Different propellants can be modeled using this formulation, with the 
difference in species factoring in to the definition of Ua, as well as in the ionization energy. The other species most 
often used in MPD thrusters are argon and lithium. Both of these propellants have much lower Alfven critical 


speeds, and much of their operating regime would then be essentially described using the Macker formula. The 
accuracy of the assumptions of constant parameters in the voltage model will have to be assessed for each of these 
propellants. Similarly, the thruster operating limit known as onset has not been incorporated into the model. Onset 
is usually associated with increased voltage fluctuations, and anode erosion. A model for this process could be 
included as a limit on u c or current in the equations. Finally, the addition of an applied magnetic field to the model 
would allow extension to both types of MPD thruster. 
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